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It is shown that the nuclear Hamiltonian can be separated rigorously as Ti = Tim + TIm- For suf- 
ficiently smooth forces the monopole part Hmis entirely responsible for Hartree-Fock selfconsistency 
and hence saturation properties. The multipole part Ti.M contains the "residual" terms - pairing, 
quadrupole, etc. - that will be analyzed in a companion paper. We include a review of the basic 
results often needed when working with multipole decompositions and average monopole properties. 
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It has not been possible, yet, to construct interactions that could satisfy simultaneously three basic conditions: 

A) to be realistic, i.e. consistent with the nucleon-nucleon (NN) phase shifts, 

B) to ensure good saturation properties, i.e. correct binding energies at the observed radii, 



' C) to provide good spectroscopy. 

As a consequence many forces have been designed for specific contexts or problems : pairing plus quadrupole , 
density dependent potentials for mean-field approaches , Landau- Migdal parametrizations for studies of the giant 
resonances direct fits to two-body matrix elements for shell- model calculations and many others. A way out 
y—{ : of this unsatisfactory state of affairs would be to exhibit an interaction consistent with conditions A,B,C above. If 
J> ' we only assume that it exists and could be reduced to an effective form smooth enough to do Hartree-Fock (HF) 
, variation, many of its properties can be discovered. The basic tool we shall need is the following separation property: 
Given a sufficiently smooth Hamiltonian ,H, it can be separated as H = Hm + Hm- Only the monopole field T-Lm 
is affected by spherical Hartree-Fock variation. Therefore it is entirely responsible for global saturation properties and 
single particle behaviour. 

I As long as we do not find ways of reconciling conditions A and B - a major problem - Tim must be treated 
• phcnomcnologically. The phenomenological enforcement of good saturation properties is quite feasible in a shell- 
model context and leads to the pleasing result that the multipole part T-Lm - which can be extracted rather uniquely 
. from the realistic interactions - has an excellent behaviour Therefore conditions A and C, as well as B and C, 

■ are mutually compatible. 

O \ An elementary argument explains the situation. The observed nuclear radii, r = 1.2yl^/'^fm, imply average inter- 
particle distances of some 2.4 fm, and therefore the nucleons "see" predominantly the medium range of the potential. 
This is a region that is well understood theoretically |^ and well described by the realistic forces. 

The picture that emerges is quite simple: the short range part of the NN interaction is not well understood, but it 
is certainly responsible for the repulsive terms necessary for saturation. Therefore, with our present knowledge of the 
nuclear forces, the phenomenological treatment of 7i„i is a necessity, but it may perhaps lead to some fresh ideas. The 
' argument is that a Fock representation that can accomodate all possible non-relativistic interactions. Although in 
principle enormously many matrix elements have to be specified, TLm is described by a small fraction of the total and 
by some formal manipulations we shall be able to isolate the very few that count (essentially the ones describing the 
bulk properties of nuclear matter) . Therefore it may possible to exhibit some simple T-Lm that contains few parameters 
and describes nuclear data satisfactorily. Since something simple in Fock space may be complicated in coordinate 
space, it may point to ways - so far overlooked - of reconciling conditions A and B. 
There are three steps in the task of giving a complete characterization of 7i. 

1. The first is technical: we must prove the separation property, and prepare the tools that make it possible to take 
the greatest advantage of the underlying symmetries: angular momentum, isospin and the unitary symmetries of the 
monopole world. 

2. As mentioned, we know much about T-Lm, but this knowledge comes in huge arrays of matrix elements. It would 
be far more useful if we could extract from this mass of numbers the truly important ones. 

3. It remains to specify the monopole field. 

We shall deal with the first point here, and with the second in a companion (next) paper |l0|| where it is shown that 
TLm is indeed simple and a hint emerges about Tim- The hint has been taken up in pd| ], a study of nuclear masses 
that provides a first approximation to point 3. 

Our main purpose now is to prove the separation property of the monopole field . The proof is far more compact 
if multipole (i.e. angular momentum) techniques are introduced. Furthermore, they are essential for the next paper. 
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French lectures [|2j are the fundamental reference on the subject and - supplemented by a good book on angular 
momentum |l3|] - they contain most of the results we may need. In assembling them in the first section and the 
appendix, we have tried to give a self-contained account, and compress in a few pages many of the things of use in 
daily shell model practice, putting some emphasis on the connection between coupled and uncoupled representations, 
and on questions of phase. 

The second section deals with separation and HF properties and it aims at making the reader familiar with basic 
technical and conceptual aspects of the monopole field. 
Notations are regrouped at the end of the appendix. 



I. BASIC COUPLINGS AND REPRESENTATIONS 



This section deals with angular momentum coupling and recoupling of operators following the techniques of French. 
There is a little here that is not adapted directly from [|l2|: some changes in notation, the use of scalar product 
whenever possible, a lesser reluctance to go over the m-scheme for simple operations and the introduction of the two 
possible analogues of hermitean conjugation in the coupled schemes. 

The reader is reminded that in the Appendix, eqs.(A16-21) contain the elementary formulae that may be needed 
in handling 3j and 6j symbols. 



A. Technical preliminaries 



The uncoupled representation or m-scheme 

We shall work in spaces containing D orbits labeled i, j, k, £ . . . each associated with operators af , a.i that act on 
a vacuum |0 > 

a+|0 >= |i > , ai\Q >^ 0,{ai,a+} ^ a,a'^ + a+a^ = (5y (1.1) 
For a given number of particles n, we can construct different states, each associated with a state operator Z~^: 
Z+{[i\,n)\0>^a+a+...a+\0>= |[i],n>, n < i2 < ■ ■ ■ < in (1.2) 

Operators of rank k take the form (we use * for Hermitean - or complex - conjugation so as to reserve the + 
superscript to creation operators) 

Rk^J2< mRk\\i']k > Z+{[i\,k)Z{[i'],k), Z([i'],fc) = (Z+([i'],fc))* (1.3) 

We shall be mainly interested in the k = 1 case, in the Hamiltonian of rank 1 + 2 (kinetic and potential energies K 
and V respectively): 

H = K + V Kij afttj - ^ Vijki afa'j'akai 

k<i 

K,,=<^K\j>; K,,,,=<^|T/|^^>, (1.4) 
and in the special fc-body operator 

z+{[i\,k)zi[i\,k) = ^ L ^ ^ ^ ( J (1.5) 

i ^ ' 

whose explicit form is dictated by its rank and the fact that it has eigenvalue 1 when acting on any state of k particles. 
The particle hole (ph) transformations 

al = di , tti = df if i <if 

al = al ^ fli = a,; ii i > if (1-6) 
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consist in interchanging creation and annihilation operators for some orbits, by interchanging vacua 

|0>^ |0>=Z+([i],/)|0>, i<if. (1.7) 
In its simplest form, f = D and all operators have to be brought to normal order. In particular 

V — ~ ^ VijMol Oj akai = — ^ Vijkeaiaja'^aj' — — ^ Vijki[a'^a'j'aiaj 

+SieSjk - SikSjg - a'^a.iSjk - a^ajSu + a^aiSji + ajajSik] (1.8) 

and 

n = D-n, J2z{[i\,k)Z+i[i],k) = J2z+{[i\,k)Zi[i\,k)=Q (1.9) 
Coupled representations 

To take advantage of the basic symmetries, operators have to be coupled to good angular momentum J and isospin 
T. The m-scheme labels i, j . . . will be replaced by pairs (rr^), (s, Sz) ■ ■ ■ where r specifies the subshell to which the 
orbit belongs, and rz its projection quantum numbers. Following French p^ , we shall introduce a product notation 
in which expressions are independent of the coupling scheme. In jt formalism, for example, a single tensorial index 
will represent pairs in spin-isospin ( JT) space 

r = jr, Tz = MTz, r = jr^, rz=mrTz^, etc. (1.10) 

(Note that r as tensor index has not exactly the same meaning as label of shell r (i.e. r = jrPr, where Pr is the 
principal quantum number). No confusion can possibly arise from this convention.) 
Expressions involving these indices will stand for products as in 

i- f' - i-Y - {-y-+'/^ [r] -2(2> + l), [F] = [JT] = (2J+1)(2T+1), (1.11) 

more generally: 

f/(r space) = U{J space)f7(r space) , (1.12) 

where U may be some 6j or Clebsh-Gordan coefficient or similar functions as in 

< 77z7'7^|rr, >=< jmfm'lJM > < ttzt'QTTz > (1.13) 

In j formalism, also called neutron-proton (np), we do not couple explicitly to good T, the tensorial indeces refer 
to a single space , and the identifications are F = J, r = jV, H = (2jr + 1) etc. (Note that when used as label, r 
must specify whether the shell is a neutron or a proton one). In is and 1st formalisms, we have {L){S) and {L){S){T) 
product spaces respectively, and the necessary identifications are straightforward as long as we do not recouple L and 
S explicitly to good J. 

Let us introduce the coupled and normalized two body state operators written in terms of ordered pairs (a^a^, i < 

■^'rr.M = C/ < rrzSSz\TTz > a+^a+^ , ^rr,(rs) = {Z+^ {rs))* 

('■.) 

a+rAs. = E < ^^-^^-|rr. > Z+^^ {rs) , Cr. = (1 + 6rsr'^^ (1.14) 

(r) 

(r^) means that the sum is restricted to < if 6rs = 1 
(F) means that F must be such that (— )'"+2'" = — 1 if (5^^ = 1 
We can always relax the ordering through 

Let us make the following identifications 

i = (r, Tz), j = (S, Sz), k EE {t, tz), I = (u, Uz)-, (1.15) 
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and take advantage of the ordering in eq.(1.14) to write directly 



t<",(r) 



(1.16a) 



(r) 



(1.16b) 



(1.16c) 



'^turs - Kstu - -{^y " ysrut — (") 



r+s-r^/ , _ ( _\r+s+t+uyT _ _( _\t+u-T^^T 



(1.16d) 



Rotational and isospin invariance are made manifest through the scalar products, which in angular momentum ma- 
nipulations must be rewritten as zero coupled pairs. There are two natural ways of doing so, depending on the order 
of the couplings: 



(1.17a) 



(1.17b) 



where Z has been replaced by the most convenient tensor, which is either Z"" = Z+, the conjugate of Z~^, or 
Z^ = Z+, the adjoint of Z~^ . The existence of these operators follows from the commutation rules of J and T with 



any irreducible tensor P^^ 



(PIT = {-V+^'Pl^^ = (_)7-7.p' 



7^ ■ 



(1.18) 



To stress the unavoidable ambiguity in sign, we have taken the unusual step of introducing two operators, the adjoint 
P and conjugate P, that are the inverse of one another: 



P7 — (— )27p7^ P7 = P7 = (_)^7p7^ p-y ^ pi ^ 



— P7 



(1.19) 



jf (_)27 _ 2 the two operators are identical. If (— )^''' = — 1, it becomes impossible to ensure - with a single definition 
- the same sign for scalar product and zero coupling of operators. Which version we choose depends on the problem 
at hand. For normal ordering as in eq.(1.17a), and in particular 



Ur — a^ ■ ar — [r]^/^(a+a~)° , 
the natural choice of basic tensors is a^ and a~ , for which we adopt the notation of French: 

With antinormal ordering as in eq.( 1.17b| ), needed in ph transforms we find 

Zr ■ z+ = [r]V2(Zj:z+)" ^ [vY'^iz+Zrf 

which simply amounts to a notation 

z^P — r ' r — r ' 

as the bar over the operators means ph transform. But it also means adjoint, and indeed: 

Zjt = by definition and Z~ = Z+ — Z'^ by definition of Z"" and by the last equality in ( 1.19| ). Therefore the 
ph transform of a state operator is its adjoint. The antitransform of a state operator is its conjugate by exactly the 



(1.20) 
(1.21) 
(1.22) 
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same argument. Note that for the uncoupled operators ph transform is Hermite an co njugation, better represented by 
its dual analogues in ( 1.18 ) than by the arbitrary choice of one of them. From ( 1.22| ) we have the basic transforms: 



The coupled operators quadratic in A and B are 

X+r^ (rs) = {ArAs)l^ , Xrr^ (rs) = (B.B,)^^ , SI (rt) = {ArBt) 
Obviously Z^-p (rs) = (rs^rr i'''^)- From the easily proved identity 

(F^f = (_)7+7'-r(Q7'p7)r 
(equally valid for conjugation), we obtain 

X+{rs) = -Xrirs) S^(ri) = {-f-^-^ Sl^^{tr) . 
For reduced matrix elements we use Racah's definition 



a 

-a. 



7 P 



< a\\P'^\\/3 > 



For any operator P'^ it is true that 

< aa,\Pl\m. >=< f3P,\{Piy\aa, >* 



and by applying (1.27) to both sides it follows that 

< a||PT||/3 (-)"-/3-7 < ;3||F^||a > , 
where we have omitted complex conjugation on the rhs because our reduced matrix elements will be real. 



(1.23) 
(1.24) 

(1.25) 
(1.26) 
(1.27) 

(1.28) 
(1.29) 



The co upled form of a rank 1 operator is deduced from the uncoupled one in (1.3) by using (1.27) and the definition 
of S"' in (|L2|): 

K^Y. < ^^^\K\*i^ > = E < ^ll^"'!!^ > i^y'^^Slirt) . (1.30) 

r,t,rz rt 

We can always rewrite an arbitrary in terms of the symmetric (S) and antisymmetric (A) operators 

Slirt) = Slirt) + i-y-'Slitr), (rt) = Sl{rt) - (-)'-*5^Jir), (1.31) 

Then calling R'^^ =< r\\K<\\t > (7)~l/^ the expansion becomes 

= ^ E [i^rt + {-r-'Rl)Sl (rt) + [RJ, - i-r-'Rl)Al (rt)] (1.32) 

r<t 

The {—Y~^ phase ensures that the spherical harmonics are symmetric if Yf*„ = {~)'^Yi^rn (Condon and Shortley's 
choice) , which leads to < r||y£||< > = (— )'^~* < > ^-^d vanishing of the A term. The phase convention 

should be changed to (-)'■- if Y£ ^ /Y«. Then F/^ = (-)^+™y£_,„, and < r||Yf||<>= (-)^-*-^ < t||y^||r >. 
Note that now Y« • = [lYl'^iYiYif instead of the usual • = {~Y[£Y/'^{YiYt)° but J • J = -V3(JJ)° and 
T -T = -V3{TT)° always, since J —> iJ makes no sense. Still, positive definite zero coupling for tensors of integer 
rank could be obtained with the change < £m£ — ra|00 >= (— )''^™[£]~^/^ to (— )'"[£]~-^/^. As we have seen, for half 
integer rank, no convention will ensure that zero coupling is always definite positive. 

Recoupling. 

Let us consider a ph transformation (1.6) and using the identification ( 1.15| ) assume that orbits j = {ssz) and 
k = (ttz) are transformed while i and £ (i.e. r and u) are left untouched. Then we have 



1/ = ^ Vijki afakO^Ui - ^ VijjiafuiSjk- 



(1.33) 



i<3 
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Now use (1.16b) , relax the ordering for projections in the contraction in (1.33) , and introduce the A B tensors, 
permuting the middle ones 



< rr^sSzl^Tz > < sSzUUzlTTz > aX^an 



(1.34) 



For clarity the internal couplings are indicated by different notations: ( )^ and [ Recoupling through a normalized 
9j symbol yields 



BtAs Bu 



nr^ ( r s T 

7 I 7 7 



(1.35) 



The contraction in (1.34) can be calculated to be 



< rrzSSz\TTz > < uUzSSz\TTz > (-)"+^-r = S^^Su^rA-) 



1/2 



(1.36) 



obtained by reorganizing the couplings with the help of 3-j symbols. All these operations can be summed up by 



-(-) 



u+t-r 



1/2 



(1.37) 



and its inverse 



+(-) 



u-t+7 



1/2 



'Jsti-'r-ii 



(1.38) 



where the first term can be written by inspection but the contraction needs some care. The hybrid procedure we 
have adopted - of mixing , a and A^ B - points to a paradox: French's notation a nd te chniqu es si mplify many 
complex coupling problems but complicate a few simple ones. (Hint: try to obtain (1.37) and (1-38) using the - 
impressive - artillery in Hsu's appendix to [|2|; a good exercise). 



B. The V and lj representation of Ti. 



Equations (1.37) and (1.38) were derived to prepare for a ph transform of some orbits. However, beyond the possi- 



bility of changing vacua - which may be quite useful occasionally - we are interested in the different representation(s) 
of the Hamiltonian that these operations entail. Accordingly, we shall always keep the untransformed notation a+, 
and Ar B^ for the orbits. 

First we consider the case in which all orbits are transformed. Introducing Sm — ci^ ■ ciu-, the coupled version of 
eq.(1.8) becomes, 



^ = - E VrstuMtu) ■ Z+{rs) - ^ im 



r 

rsrs 



(1.39) 



r<3 

t<u,r 



r<s,r 



r 

rstu 



[r] [s] J [r] [s] 
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All contractions are calculated using ( 1.36| ) 



Next we examine the transformations associated to ( 1.37 ) and ( 1.38 ) in which only the middle operators are 
interchanged. It is convenient to allow for complete flexibility in the summations and we shall relax the restrictions 
T l£ s, t < u hy replacing the factors by the P^s convention 



(1 + SrsV^'^ 

Prs = (1 + Srs)^^^"^ ii r < s , Prs = ^- if uo restriction 

so that the sums could be interpreted as restricted or not restricted. We write therefore 

'■<!> (rstu)r 
t<u,r ^ ' 

and from now o n we set Zyijs) = etc. 
According to (1.37), V can be transformed into 



(1.40) 
(1.41) 



where 



(r) 



ut ^ 



u t -f 



(Remember that X](r) means that we sum over Pauli allowed T). Eq.( 1.3§| ) suggests an alternative to ( 1.42 ) 



{rstu)-y 



1/2 



(1.42) 

(1.43a) 
(1.43b) 

(1.44) 



where each multipolarity 7 is associated with a two body operator. The obvious check that (1.44) is indeed ( 1.42| ) 
comes from 



\f + 1 — s 



1/2 



— [s] ^ ^russ ■ 



(1.45) 



The proof is left as an exercise (use ( 1.43a ) and Racah sum rule (A. 19)). 

In a jt representation, by introducing explicitly the isospin in F = JT, 7 = At, r = etc and the 6-j values 

ri/9i/9T^i 00 01 10 11 

1/2 1/2 r = -1/2 1/2 1/6 



we find 



AO 



,A1 



iJ) 



■E(-)^' 



iJ) 



ju jt 



and reciprocally 



V, 



JO 



(1.46) 

(1.47a) 
(1.47b) 

(1.48a) 



71 



4 EH 



(1.48b) 



Wh en the Hamiltonian is written as in (1.41) we speak of the normal or T^-representation while we call the form 
(1.42) multipole, or w-representation. 
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II. THE MONOPOLE FIELD TLm AND THE SEPARATION n = nM+ Hm 



Before going into the technical problems, it is worth explaining how the monopole field appears. Several of the 
statements that follow need a formal proof. It will be supplied in the body of the section. 

The separation of the Hamiltonian into an "unperturbed" and a "residual" part , TL = TL^ + 7i', is at the heart of 
many-body physics, and the idea that Ho, could be represented by some central - i.e. single particle - field is of great 
heuristic and qualitative value. However, to decide whether it makes sense quantitatively, we must understand how 
7Y, of rank-2, can be approximated by a rank-1 operator. 

What can be done cleanly is to define Tio as K plus a two body part so that 

n^K + Y^Vrs nrins - 6rs)/il + 6rs) + H' . (2.1) 

Tio has two equivalent properties: 

Complete Extraction. When written in a multipole representation 7i' contains no number operators. It means that 
all the A = terms that are diagonal in the basis we are using have migrated to Ho- 

Trace Equivalence. The expectation value of Tip for any basic state is the average energy (i.e. the trace) of the 
configuration to which it belongs. (A configuration is the set of states with fixed number of particles in each orbit). 

To recover and examine the notion of central field we choose as vacuum some determinantal state, i.e. we set 
Tir — Dr bel ow the Fermi orbit and Ur — Q above. The single - particle and single - hole energies are then calculated 
from cq.( |2.l|) , by adding or removing a particle from the vacuum. In a HF calculation these energies play an important 
role and invite the interpretation that the determinantal state generates a mean field in which the particles move. 
This interpretation is very good when there is a dominant agent in the Hamiltonian that is responsible for most of 
the field, as in atoms and planets. 

In nuclei it is simply wrong. What the HF calculation produces is a basis of orbitals. When Ti is written in 
that basis, Tip can be extracted, and it will yield indeed the HF value for the vacuum and single fermion energies. 
However, there is no reason to linearize Ho (i.e. to approximate it by a central field) when estimating energies of 
other configurations: the quadratic effects grow fast. To fix ideas: typically Vrs ~ 300 keV in the pf shell, not a 
large number, but it multiplies the UrUs operator that may become 0(10) or even O(IGG), and drastically modify the 
effects of the central field. 

There is much empirical evidence that at fixed total number of particles it is a good approximation to keep the 
same basis for all orbits in the vicinity of the Fermi level, and therefore Hq is fixed. However we also know that, when 
we change the number of particles, the orbitals and therefore Tio must evolve. 

Since the evolution of the orbitals is associated to unitary transformations of the basis, the rir operators become 
linear combinations of S'J'^j ones. Therefore we must generalize the definition of TIq, by extending the notion of 
"complete extraction" to all A = operators, not only the diagonal ones. The resulting object is what we call Hm- 
As the separation Ti. — Tim + Ti-M is invariant (representation- independent), Tim provides the - mathematically and 
physically - natural definition of unperturbed Hamiltonian. 



A. Separation of Tim 

Although there is only one Hm, it will have different aspects in jt and j formalisms. We differenciate them by a 
tag: 



Hr. 



Tin 



TirnT is constucted by extracting all the possible 7 = 00 and 01 contributions to eqs.( 1.42 ) or ( 1.44 ), while Hmnp 
contains all the possible 7 = terms. O bviou sly the kinetic energy is part of Hm, whose rank 2 component must have 
the form of the A = contributions to ( 1.44 ) which we call V\=q. To give them a transparent aspect we replace S'™ 
and S^l by 



(2.2) 



which, for Sr 



1, reduce to Sr 



Tr- The two body operators 



Srtsu — CrsCtui^SrtSsu ^st.Sru^ 



(2.3a) 
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^rtsu — CrsCtui^'^rt ' '^su stSru) 



(2.3b) 



in turn become for 5rt — 5su = 1, 



nr{n - (5rs)/(l + Srs) , 



(2.4a) 



^rrss — (-^r ' ^^rf^rs)/(l f^rs) • 

and the A = contribution to V can be written as 



(2.4b) 



(2.5) 



which follows directly from ( ^.3a ), (2.3b) and (1.44) provided we remember that Tj-t ■ Tsu = — \/3(TrtTs„)°. 
Although the full ra nk-2 c ontrib ution t o HmT has the form Va=o, the correct values of the ujrtsu ''^^'^ '^rtsu cannot 
be extracted from eqs.( 1.47a ) and ( 1.47b ) as 



(J) 



(2.6a) 



,01 



-1/2 



(2.6b) 



There are two reasons. One is the Pauli principle, which operates when r = s or t — u, i.e. when 6(P) — 1 with S{P) 
defined as 



S{P) = 1 - {1 - 5rs){l - 5tn) 



(2.7) 



A direct indication of the problem comes if we consider a single shell with a constant interaction, T^^^^,^ — 1. Then 
from (1.5) 



Z+p • ZrrT = ririUr - l)/2; 

but from ( ^.6a| ) we obtain woo — {Dr — l)/2 with Dr = [r], and inserting in (2.5), 

^A^O = {{Dr - \)l1Dr)nr{nr - l)/2 

, indicating that we miss the correct result by over a factor 2. The remaining half must be found among the terms 
of Vx/o J which in the case of (5(P) = 1 cannot be linearly independent since they are twice as numerous as the 
ones. 

The other problem is related to exchange terms that show for jV = js = jt — ju and 5{P) = 0, i.e. when 5{e) = 1 
with 



6{e) = {5rs - Srs)iStu - Stu) , 



(2.8) 



Then, setting = 1 in (2.6a) and ( 2.6b| ) and remembering that V^^^f = — (^ )'^'''^Kstti (since (— )^'' = 1 in jt 

formalism) 



-[rs]-'^'J2^J]{~y{KiL--3Kil) 
{J) 



(2.9a) 



(2.9b) 



indicating that much A = monopole strength is hidden in the i^rLu terms with A 7^ 0. 
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1. The form of TL„ 



To achieve complete extraction we specify Ti-mT by defect^ as suggested at the beginning of the section. It is Hm 
that will be defined so as to contain no Va=o terms. To illustrate how the prescription works we consider first the 
(5(e) = case which is simp ler. 

Since TLm has the form (1.41), let us replace the V^^^ matrix elements by 



* ' rstu 



y/sL - yrstJrtSsu if ^(e) = 0, 



with the "centroids" defined so as to make wJ;J^„ vanish: 

(J) (J) (J) 

Then TimT is automatically defined in the V representation as 



(2.10) 



(2.11) 



n,nT{S{e) =Q)= K+ SrtSsuV^stu ^rsJT ' ^tuJT 

= . (J) 



(2.12) 

The remaining task is to write the operators in TimT in terms of Srtsu and Trtsu, which we do easily by calling 



upon eqs. (1.48a) and (1.48b), where {5rs5tu = 1 of course) 



which means 



1-1/2 



-1/2 



(2.13a) 



and inverting 



J 



tujl 



(2.13b) 



(2.14a) 



(2.14b) 



The exchange contribution to TimT is a bit trickier to extract, because when (5(e) = 1, two extra operators have to 
be considered. They can be calculated from (2.14a) and (2.14b) by exchanging t and u, 



(2.15a) 



and by combining (|2.14a|) , ( |2.14b|) , ( ^.15a| ) and (^.15b| ) 

(1± (-)•') _ 1 



(2.15b) 



(2.16a) 
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rtsu 



(2.16b) 



Now it is possible to specify completely the matrix elements of T^m, (^(e) defined in (p.^) 



(1 - 5(e))l/^^t„ + 5{e)V^,t^ , p = sign (-) 



rstu ~ ''rstu "rfJsu 

from which it follows that Ti-mT rnust be 



(2.17) 



rstu 



f <Ti,T,p=± 



(1±(-)^) 



(2.18) 



Through eqs. ( p4a| ), ( pl4b| ), ( p6a| ), and ( pl6b| ) we can obtain the form of Ti.„iT in terms of the monopole operators 
by regrouping the coefficients affecting each of them. To simplify the presentation we adopt the following convention 

a = rstu r < s,t < u, Sj-tSsu = 1, BUT 



then 



H™t + ^(1 - 6{e)){a^Sa + 

a 

+Sie)iaiSa, + biT^ + a^Sa, + b^T^), with 







a-a = 




+ vS) 






ba 






»i = 




+ 


+ + 




"i = 




- 


-v+° + 

^ a ' 




K = 




+ 


^ a 


^ a } 


K = 




- 

*^ a 


' ^ a 


^ a. / 



{J) 



From ( 2.11 ) and its generalization to (5(e) = 1 we find for the ccntroids 

{J) {J) ^ ' 

vtL = E y/sUm ± (-)'')/Ew(i ± (-)') 

J .7 
J 

Dr = H, 5rt5su = 1 always; (5(e) = 1, for V^^^ . 



(2.19a) 



(2.19b) 



(2.19c) 



(2.20) 



2. Form of Hmnp and Hmo 

In j formalism Jimnp is 'HmT under another guise: neutron and proton shells arc differentiated and the operators 
Trs and Srs are written in terms of 4 scalars 5^^^ ;x,y ~ n or p. 
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We may also be interested in remaining in jt formalism and extract only the purely isoscalar contribution to TimT, 
which we callTimO- The power of French's product notation becomes particularly evident here, because the form of 
both terms is identical. It is left as an exercise to find 

n,nnp or H„,o = K + Y, [vo.Sc.il - (5(e)) + ]^{{V+ + l/-)5„ + [V' - V+)Ss.)5{e)\ (2.21) 

with 

(r) (r) 

Y,[T] = Dr{D,-5{P))l{l + 5{P)) 
(r) 

r 

^[r](l ± = DriDr T i-ri (2.22) 

r 

Of course we must remember that for Hmnp Dr — 2jr + 1, (— )^'' = —1, F = J, etc., while for Timo 
Dr = 2{2jr + 1), (-)2'' = +1, F = JT; etc 



It should be noted that Hmo is not obtained by simply discarding the b coefficients in eqs.( 2.19 ) , because we can 
extract some 7 =00 contribution from the Ta operators. The point will become quite clear when considering the 
diagonal contributions. 

B. Diagonal forms of 7i„i 

We are going to specialize to the diagonal terms of HmT and Hmnp, which involve only V^^^^ matrix elements whose 
centroids will be called simply Vrs and V^J^ (the overline in Vr st,,, w as me ant t o avoid confusion with possible matrix 
elements T^rstu oi' Kstm it can be safely dropped now). Then ( ^.19 ) and ( ^.21 ) become 



^mT =K + Y^i^rs nr{ns - drs)/{l + Srs) + brsiTr ■ - |n(5,,)/(l + 5rs)) (2.23a) 



r<s 



^tnp or nio =K + J2yrs nr{n, - Srs)/il + Srs) (2.23b) 

r<s 

We rewrite the relevant centroids incorporating explicitly the Pauli restrictions 

Dr{Ds-Srs) " DriDs+26rsi-V) ^' ' 

^rs = \ (3Ks + ) ^Vrs + \ Ks = ~ (2.24b) 



The relationship between ars and Vrs (which should be checked) makes it possible to combine eqs.( 2.23a ) and ( 2.23b| ) 
in a single form 

nt^K + Y.^-- "'■("« - ^-)/(i + + ■ - 4(j^''-i) '^^^j + ^'■^^ (2-2^) 



in which now the hrs term can be dropped to obtain TYJ^q or Tifnnp- Note that the diagonal terms depend on the 
representation: 

T^mnp 7^ i'^ general. It is useful to have at hand the relationship between matrix elements and centroids in 

the jt and the j schemes. A single shell r in the former becomes a pair r„ and Vp in the latter 
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*'rnSptnUp y^rstu ^ ^rstuJ 



(l+(5,,)V2(l + 5,Jl/2 



2 



T/ T/'^l 

"rpSptpUp — <'rstu 



and for the centroids {x, y — n or p, x ^ y) 



(2.26) 



(2.27) 



1. Trace equivalence and propagation 



We have already mentioned the important - and known p^- property of that its expectation value for a given 
state reproduces the trace of Ti. over configurations Q = {ni n2 . . .rir . . .) or Q = {niTi 712X2 . . . rirTr), to which the 
state belongs. The notation Q stands for fixed or UrTr in each subshell. Averages at fixed seniority and for 
other symmetries are also of much interest [1l6| . Here we propose a very general proof that relies on Schur's lemma 
and on a trick used by French ||T^. Let now Q stand for any symmetry associated with states \YQ > oi n particles 
(n is included in Q). The trick is to construct \YQ > by acting on |0 > instead of |0 >. Let us introduce Vq of rank 
n = D — n, 



T^Q^H^YQ^yQ^ \yQ>=ZYQ\^>, \YQ>^Z+^\Q>; 



(2.28) 



The quantum numbers Q are associated to the ph transforms of the operators e.g. Q — {niTin2T2 . . .). Vq is a 
projector 



Vq\YQ' >^ Sqq,\YQ > if n'< 



(2.29) 



From Schur's lemma it follows that it can be written in terms of the Casimir operators associated with the Q 
quantum numbers. It is seen that Vq is a generalization of Vk in (1.5) and for Q — {ni 712 . . .rir ■ . ■), it is obvious 
that 



Vr 



n 



ni 712 ■■■nj.. 



n 



ni 712 ■■■'rir 



(2.30) 



where we are forced to distinguish the operator hr and its eigenvalue n^. Note that Vq is not the ph transform of 
Vq. 

We call dg =J2y YQ\YQ >~< 0|I?q|0 >; the dimension of the space over which we take the trace < H >q; 
and WxxqZ^qZxq the diagonal contributions to TL (in our case x = J, q ^ rsT). Now we can calculate the trace (a 
step by step explanation follows the equations): 



< n >Q- d-Q^ E < E WxxqZtqZxq\YQ >= 



(2.31a) 



dg' J2 < 0\Z+^Z+ZxqZyQ\0 >-- 

xqY 



(2.31b) 



(2.31c) 



(2.31d) 
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(2.31e) 



(2.31f) 



Line by line we have 

a. definition of trace 

b. insert the definition of \YQ >= ZyQ 



|0> 



c. French's trick of using Zyq\0 > instead of ZyglO > for \YQ > allows to permute operators. Then use eq.( 2.2S| ) 

d. the crucial step: since the expectation value of Vq is independent of x wc take it out of the sum for some 
arbitrary x' . 

e. by the same token we sum over x' and divide by the number of terms, then commute again the operators in q 
and Q. 

f. rearrange e to make clear what the trace equivalent operator is [p^ . 

If we specialize this result to our problem, eqs.(2.14a) and (2.14b) for r = i, s — u provide Vq and eq.( ^.ll ) provides 
J2 Wxxq- Therefore, Ti.'l^rp is indeed the trace equivalent of H. Reference to Schur's lemma may be avoided for 

our particular problems. What has to be show n is that we have exactly the number of operators required to calculate 
Vq. For 7i^„p the problem is solved in ( 2.30| ). For Ti^T ^^^^ form of Vq is not so simple, but we can easily check 
that the number of operators available meets exactly our needs. Let us do it for 1 shell (or for total n and T), so that 

Va=Vn 



■■q — i^nT- 

For n = 2 we have from eqs.( ^.14aD and ( |2.14b| ) , (T^ = T{T + 1)), 



Vo 



1 



{n{n 



V21 



1 



(3n(n - 2) + AT^) 



(2.32) 



For n — 3 we need V^i and V^^ which we construct out of n'^^^' and (n — 2) (T^ — 3/4n), the two available rank 3 

operators. For n = 4 we need V4Q, V41, V42 and we have n^, n^T^T^ from which we can extract rank 4 operators 
and work out 



P40 = ^(T^ 



n{n + 2) 



2^41 = 7^ T 



n{n + 2) 



Vi2 



T'-]n{n + 2) 



-(n-2)(n-4) 



24 



V 



40 



V 



41 



(2.33) 



Obviously the construction can be extended to any VnT 



If we want to deal with averages at fixed intermediate T values e.g. 

Q = {niTi n2T2 nsT^ T12 T), even the counting may become hard: Schur's lemma is not trivial. Still we can do 
without it, by direct averaging through standard multipole techniques, as in 



J2 < ?^l^lyln2^2r2^|| [{AiBAr{A2B2r]° \\n1r1Y1n2r2Y2r >^ 
Y1Y2J1J2J 

'riFaF' 



E 

YlY2.hJ2J 



(-) 



ri+r2-r-7 



raFiF 



< riyi|i57iiriyi >< FaFsll^JlirsFs > 



(2.34) 



where, using the sum rule (A19) for J, only 7 = 01 and 00 survive, which are absent from Hm, establishing Hm as 
trace equivalent for these operators. The only averages that cannot be done in this way are the basic ones involving 
n and T only, but these we have already. 

The idea to associate a trace to an equivalent operator and not only to a number, leads to French's idea of 
propagation. It stresses the fact that by combining the general properties and particular symmetries of a given trace 
equivalent operator we may determine it at some convenient and limited set of Q values and then propagate i.e. apply 
it to any Q. The task is hard in general but for purely isocalar averages it is often possible to to guess outright the 
solution, as for the "width" of Hm- 
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a\n)^<nl,>n=^^^^^a'{2), a^^) = ^Y.Hm (2-35) 

since the trace equivalent operator must be of rank 4, reproduce cr^(2) at n = 2, and vanish at n = 0, 1, n = 0, 1 
because the energies of the corresponding states are entirely given by Tim ■ 



C. The Hartree-Fock property of Tim- Strict rank operators 



Equation (1.39) contains all the possible contractions that may appear in any possible change of vacuum. By 
construction they all vanish for T-Lm since their coefRcients are proportional to centroids that have all migrated to 
Tim- As a consequence Hm is of rank 2 both for particles and holes (we call this property, strict rank-2) and its 
expectation value vanishes for any closed shell and any single particle or single hole state, a set we denote by "CSib 1" . 
As a consequence, 

<cs± i\n\cs ± 1 >=< cs ± i\n,n\cs ± i > (2.36) 

which remains valid under unitary transformations ai^ = Urxb^ because centroids become linear combinations of 
centroids. Therefore, if they all vanish in one representation, so they do in any other representation and Tim remains 



a strict rank-2 operator. Of course at any particular basis we can replace 7i„i by in (2.3(:) 
From these remarks it follows that 

Hartree-Fock variation ofTi. is Hartree-Fock variation ofTi.- 



Eq. dlH) can be obtained more directly by arguing that any CS ±1 state is the only member of the configuration 



to which it belongs. If we have invoked eq.(1.39) instead, it is to introduce the concept of strict rank-2 operators. 
It can be certainly generalized and it is quite useful for monopole operators, although it must take a weaker form 
since for Hm the property guarantees v anish ing for all CS ± 1 states, which cannot be the case for any diagonal 



contribution Tim . For instance: in eq. ( 2.25 ) the isoscalar part has been cleanly separated because T — 3/ An of 



rank 2 in eq. (2.23a) has been changed into — 4( ^ "j^-, which is of strict rank 2 for doubly closed shells. 



For number operators it is easy to check that 



^Dr DiJ \Dr Di 

is of strict rank 1 for closed shells containing 1 and r, while 

'o (2) (2) \ n n 

ZUrU^ ni. n\ \ JJrL)^ 



, nr rii \ I nr Ut , , ^ 

^r= J^-^]Dr = - -^--^ Dr = -Tr (2.37) 



Trs I 7-» (o^ ro\ 1 r. — ^rs — ^rsi'^ri '^s) (2.38) 

is of strict rank 2 for closed shells containing r and s. 

It is seen that the set of operators ni, n2, . . .n^ can be replaced by the total number n and v — 1 operators F^ (in 
which the choice of ni is arbitrary). Similarly we have n'^', nVr and F^s with r 7^ 1 and r ^ s, replacing the nrUs 
set. 

The form of TYJ^jq obtained with this construction is given in where it made possible an efficient monopole 
phenomenology (beware of two erratae in this reference: in the definition of read Nr not Nir, and in that of Wr read 
Ns not Nr). For the form of Ti-mnp refer to [|l^, where several applications are given. In particular a reinterpretation 
of BCS as a number conserving theory which hinges on what is the fundamental point of the construction: That it 
singles out n - always a conserved quantity, and the only one of strict rank - from all the other number operators. 

To understand the interest of this result in giving a full characterization of TLmo we note that the energy of a closed 
shell to which the nF^ and F^s operators do not contribute, is given by Wn{n — l)/2, with 

W = Drs = Dr{D s - 5rs) / {I + 5rs) ■ 

It means that a single parameter W is in principle sufficient to describe the energies of all closed shells. The catch 
is that all the strict rank properties apply to the particular closed shell being considered and W is necessarily 
space-dependent, since defined only for a given set of occuppied orbits. In a way we have gained a lot by isolating 
the strongest contributor but it is of no use to have a parameter that is space dependent if we do not know the 



dependence. In |10|, we shall see how the problem is naturally solved, once we consult the realistic interactions. 

The separation property of Tim had been anticipated |l^, used 0, and very sketchily proved for Hmnp 0]- The 
full derivation was long overdue and it might have been given in a briefer form but we hope that the longer one we 



have chosen will be more useful. 
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APPENDIX A: 



Reduced matrix elements for ^, a, rY\ and q operators 

Conventions: radial wavefunctions are positive near the origin, 
j = + s (not s 

Condon and Shortley phases for Yum- 



A) a = 2's and i operators. 



< m\f >=< m\j >= 



j' = 1 + 1/2 J = £-1/2 /(j)-(j(i + i)(2i + i))'/' 
< >= /(/)// < jMj >= -.f{j)/{j + 1) 

< fm\f >= u{f)/f < m\j >= (£ + + 1) 

2£(£+l)"l^/^ 



2^+1 



= -<m\j>=-<mf> 



B) operators. ^ = 



2i+l\ 1/2 



( 3 3 

<pmr^cx\\p'ej' >=<pt\r^\p'e > {-y'-^-y^\jj']y^ 

Vl/2 -1/2 



Bl) r^Yi operators. 



<pe\r\p + ii + i>= \i-{p + e + 3) 



< pe\r\p + 1 e - 1 >= -\i -{p - e + 2) 



<pij\\rCi\\p + i£ + ij + l>-- 
< pej\\rCi\\p + 1 e + 1 j >= 
<pej\\rCi\\p + ie-ij > 

<pij\\rCi\\p+l£-lj~l>= 



p + £ + 3 



(2j + l)(2j + 3) 



p + £ + 3 



p-i + 2 



2{2j + 2) 

2j + l 



2j-(2j+2), 
2j + l 



2j(2j + 2) 



1/2 



1/2 



1/2 



p-i + 2 



(2j-l)(2j + l) 



4j 



1/2 



B2) r'^Y2 operators {qohu only) 



<pe\r^\pe >=p + 3/2 <pe\r'^\pe + 2>= -{ip-£){p + £ + 3)y^^ 



<3nr^C2\[j£>-- 



p + S/2 



(2i + l)(2i-l)(2i + 3) 



<j£\\r''C2\\j + l£>=ip + S/2) 



2i(2i + 2) 

3 (2j + l)(2i + 3) 



1/2 



1/2 



<j£\\r''C2\\j + l£ + 2> 



2 2/(2.7 + 2) (2i + 4) _ 
3{p- £){p + £ + 3)(2j + l)(2j + 3)T 



<i^l|r'C^2||j + 2£ + 2>=- 



2 2j(2i + 2)(2i + 4) 

3 (p - + £ + 3)(2j + 5)(2j + 3)(2j + 1)^ 



(2i + 4)(2i + 2) 
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Useful formulae involving 3j and 6j symbols 

= {^y+i'-J < fm'jm\JM >= {-y+i'-J < j - mj' - m'| J - M > 



(A16) 



TO to' -My/ ^ l^m' TO -My ^ l^-m -to' My ^ ' 

The 3j symbol is invariant under cyclical permutations. 

< jm j - to|00 (A18) 

EH'w{/ ^' fc} = (-F+^'[j/]'/''5..o (A19) 

fc} = (-)^+^'+'^by]-^/^ (A21) 

The 6j symbol is invariant under permutations of columns and under interchange of upper and lower indices in two 
columns. 

NOTATIONS 

We collect here a list of notations used throughout the text, indicating, when useful, the equation or section in 
which they are introduced or explained. Occasionally some comments are made. 

WARNING: The same notation may mean different things, but once this is kept in mind it should be very easy to 
give the correct interpretation in context . 
Three general points 

1) = F{F — 1) • s{F — n + 1), integer n. Boole's factorial power notation. 

2) Hardy's Order notation f{x) — 0{x) means f{x)/x < constant, for x tending to some limit (x — s- cx) or x — s- 
usually). 

3) = TkK, (never T^'^) is a tensor operator of rank k and projection k. Applies to coupled operators as in 
{T^S^ )^ , in which case the rank K is always a superindex. Whenever possible, projections are omitted. Zero 
coupling: {T^S^Y', when working in product space (JT, LS or LST), means full scalar. E.g. in JT space = 00, 
i.e. scalar, isoscalar. 

Symbols 

• =^ implies; transforms, tends to, becomes. 

• [j] = 2j + 1, but see ( |l . 1 1| ) for product spaces. Also [i] = ii, ?2, ■ • ■ in, but only in the very first lines of section I. 



< >Q = trace ( |2.3l| ) 



• • scalar product T • (S"")* = J2k T^i^^T 

• * Hermitean conjugation for operators, complex conjugation for numbers. 
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Creation symbol as in a"*", Z"'", which are creation operators. 



• — ph tra nsfor m (1.6) 
Adjoint ( PI) 
Centroid (poj) 



• ^n, J, T are used as operators and eigenvalues. When necessary to distinguish them we use n, J and T for the 
operators. 



J2( ) restricted sum ( 2.16a ) and ( 2.16b ) 



Acronyms 

CM: center of mass. 

CS, CS±1: Closed shell, CS plus 1 particle and 1 hole states. 

HF: Hartree Fock. 

HO: Harmonic oscillator. 

NN: Nucleon-Nucleon. 

np or vt:: neutron, proton. 

ph: particle hole. 

rhs and Ihs: right hand side, left hand side. 
SM: SheU Model. 

Operators, variables 

a+, Qr (1.1), Bn (1.23), creation, annihilation operators 



A: total number of particles. AJ^: antisymmetric multipole (1.31) 
d: dimension of space. 

D — J2r B'r = J2M- Degeneracies of space and orbit r. 

G: pairing strength. 

Ti.: Hamiltonian in full space. 

H: Hamiltonian in subspace. 

i, j, k, £: individual orbits . 

J, jr-: angular momentum (total or of orbit r). 

K: kinetic energy. 

L, ir'. orbital angular momentum (total, or of orbit r). 
M: nucleon mass. 

n, rir'. number of particles (total or for shell r), 
N: number of neutrons 

p, Pr- principal quantum number pr — 2i>r + ^r- 
p: momentum ( see u nder r). 
Prs convention ( 1.4C| ). 

r,s,t,u: subshells. r = (jrP r)- 

51 = S2 (rs) = (ArBs):^ (1l24|) SZ (fM - 

r,,tr7pi),(pb). 



n, rir'- number of holes. 



^rstu : 



T: isospin. 

V: potential energy. V representation section LB 
Vijki, two body matrix elements (1.16). 
H) Vrs (HH), centroids. 
W representation (1.40). 
(II.8). 
(rs) = {ArAs)l 



rstu f 

Vrstu ( |2.20D (| 



XrvSrs) = {BrBsYr^ (1.24) 



1.14). 



Ztt, {rs) or Z+p 
Z: number of protons. 
a = rstu (i.e. W^.^^ = W^) ^ ^ ^ 

7 = At, r = JT product tensor indices in the V and omega representations. 
Srs Kronecker S. 

Srs = Sjrjs SlTrT^s , T^r = (~)^''- 
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Crs = {l + Srs)-'/^ 

A, r: angular momentum and isospin in uj representation. 
Vr = nodal quantum number: pr ^ 2vr + ^r- 
V, tt: neutron, proton. 
Tr'. isos pin projection fo r orbit r. 



is pm pro.icction iq r 
(|1.42| , |l.434l-43b|) . 



(jj = HO frequency. 
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